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We report a self-consistent microscopic theory of characteristic length scales, vortex structure and 
type-1.5 superconducting state in two-band systems using two-band Eilenberger formalism. 
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I. INTRODUCTION 

The usual classification of superconductors charac- 
terizes materials by the Ginzburg-Landau parameter k 
(which is the ratio of the characteristic length scale of 
the order parameter variation ^ and the magnetic field 
penetration length A)^. The remarkable property is that 
within the GL theory of single-component superconduc- 
tivity ti determines the major features of the phase di- 
agram of the system in magnetic field. In type-I super- 
conductor K < (i.e. order parameter is the slowest 
varying field), vortex excitations have attractive inter- 
action and are thcrmodynamically unstable in applied 
magnetic field. Thus in an applied field a type-I sys- 
tem forms macroscopically large normal domains^. For 
At > 1/^/2 (type-II superconductivity) vortices are thcr- 
modynamically stable and interact repulsively yielding a 
new phase in strong magnetic fields: a lattice of quan- 
tized vortices^i^. In the Bogomolnyi limit (k = l/-\/2) 
the vortices do not interact in the Ginzburg-Landau the- 
ory. However indeed it should be remarked that going 
to a deeper microscopic level there are always "next- 
to-leading order" microscopic corrections. These cor- 
rections, though unimportant even slightly away from 
this limit, provide weak non-universal intervortex inter- 
actions when K is very close to 1/a/2 see e.g.— i^. Ap- 
parently a counterpart of this limit is also possible in 
multi-component systems. However in this case the 
Bogomolnyi limit could appear only via quite extreme 
fine-tuning of parameters and therefore is not of much 
physical relevance. In this work wc arc interested only in 
the entirely different physics of intervortex interactions 
and magnetic response of multicomponent systems orig- 
inating from the different funamental length scales very 
far from any counterparts of Bogomolnyi limit. 

A question which attracted much attention recently 
is whether the type-I/type-II classification is sufficient 
for characterizing the rapidly growing family of multi- 
component systems of physical interest^. A clear cut 
example of the system where type-I/type-II dichotomy 
does not hold is the projected coexistent electronic and 
protonic (or deuteronic) superconductivity^ in hydrogen 
isotopes, their mixtures and hydrogen rich alloys at ul- 
trahigh compression as well as the coexisting protonic 
and I]~-hyperonic superconductivity in neutron stars. 



These systems have C/(l) x L/(l) or higher symmetries and 
thus several fundamental length scales associated with in- 
dependently conserved fields. Consequently the system 
cannot be characterized by a single dimensionless param- 
eter K. In an applied field the only thcrmodynamically 
stable vortex solutions are "composite" vortices where 
both condensates have 2tt phase windings. Consequently 
such vortices have cores in both components^i^. Impor- 
tantly it also acquires a new regime^ for which the term 
"type-1.5" was coined recently^. In that regime like in a 
type-I case the characteristic core size of one of the com- 
ponents is larger than the flux carrying area. The overlap 
of these cores produces attractive intervortex interaction. 
However, in contrast to type-I case, these vortices have 
repulsive interaction at short rangesi^^"— . This kind of 
non-monotonic vortex interaction results in the appear- 
ance of the additional "semi-Meissner" phase in low mag- 
netic fields. In that phase vortices form clusters where 
because of overlap of cores the slowest varying density 
component is suppressed. Moreover these vortex clus- 
ters coexist with the domains of two-component Meissner 
state. 



The recent experimental works proposed that two- 
bandit electronic material MgB2 belongs to the type- 
1.5 case^ii^. The principal difference with the discussed 
above U{\) x U{\) theory is that interband coupling 
breaks the symmetry down to U{\) (for a recent discus- 
sion of microscopic details see e.g] ^^'^^ ). Therefore there 
is a single superconducting phase transition at a single Tc. 
However, at the same time the system has two gaps and 
two superfluid densities, which, in general vary at distinct 
characteristic length scales at any finite distance from 
Tc- Therefore the type-1.5 magnetic response can arise 
even infinitesimally far away from Tc from the interplay 
of two density modes which originate from the underly- 
ing two-gap physics. This behaviour was demonstrated in 
the framework of phenomenological two-component U{\) 
GL modelsiS— . Here we develop a theory of type-1.5 su- 
perconductivity based on a microscopic theory without 
involving a GL expansion. That is, in this work we use 
the Eilenberger formalism and demonstrate the existence 
as well as describe basic properties of type- 1.5 supercon- 
ductivity in multiband materials. 
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II. MICROSCOPIC DESCRIPTION OF VORTEX 
STATE IN MULTIBAND SUPERCONDUCTOR 

A. Eilenberger formalism 

We consider a superconductor with two overlapping 
bands at the Fermi leveli^. The corresponding two 
sheets of the Fermi surface are assumed to be cyhndri- 
caL Within quasiclassical approximation the band pa- 
rameters characterizing the two different sheets of the 
Fermi surface are the Fermi velocities Vpj and the par- 
tial densities of states (DOS) I'j, labelled by the band 
index j = 1,2. We normalize the energies to the critical 
temperature Tc and length to ro ~ hVpi/Tc- The system 
of Eilenberger equations for two bands is 

VFjnp (V + iA) fj + 2a;„/j- - 2Aj.gj = 0, (1) 
VFjn^ (V - lA) /+ - 2uj,J+ + 2A*gj = 0. 

Here w„ = (2?i + 1)t:T are Matsubara frequencies and 
vpj = Vfj/Vfi- The vector Hp = (cos6'p, sin^p) parame- 
terizes the position on 2D cylindrical Fermi surfaces. The 
quasiclassical Green's functions in each band obey nor- 
malization condition g'j + fjfj^ = 1- The self-consistency 
equation for the gaps is 

X^tJ2 KfjdOp. (2) 

The coupling matrix \ij satisfies the symmetry relations 
niXi2 = n2\2i where rii are the partial DOS normalized 
so that ni + n2 = 1. We consider An > A22 and there- 
fore refer to the first band as "strong" and to the second 
as "weak". The vector potential satisfies the Maxwell 
equation 

VxVxA=j (3) 

where the current is 

j = -T ^ aj^Im Upg^dOp. (4) 

j=l,2 n=0 

The parameters ctj are given by 

CTj = TT ( — ) {rQVFl)'^l^jVFj- 



B. Multiple masses of the A fields 

First we focus on the structure of an isolated axi- 
ally symmetric vortex characterized by the non-trivial 
phase winding of the gap functions Ai^2 = |Ai^2|(?')e"^. 
We begin by finding the asymptotics of the gap func- 
tion modules |Ai.2|(f) at distances far from the vor- 
tex core. In this case the Eilenberger Eqs.([T|) can be 



linearized by generalizing the methods used for single 
band superconductorsi^. The details of the asymptotics 
derivation arc given in the AppendijAl We rewrite the 
Eqs.([T]) in terms of the deviations from the vacuum state 
values A J = Ajo -\Aj \ and fj = fjo - fj, /+ = f^^ - // 
keeping on the left side the first order terms. Then we 
take the real part of the EciS.(IT]) to obtain the following 
system 

VF,r,pVr^^ + 2u;„/;. - X£^- (5) 
VF.npVf,^ + 2^/^^. - ^A, = X:^^, 

where 17„, = ^cj^ + A^^, f^^ = Re[f, + /+] and = 

Re[fj — /j^]. In Eqs.® the higher order terms in A^-, 
/ and /"'' are incorporated in the right hand side (r.h.s) 
source functions ^E(d)i = -'^E(d)j ("p, w„, r). 

The solution of Eqs.(l5]) can be found in the momentum 
representation /|,^^j^.(k) = / /^^^^^.(r) exp(-ikr)d2r. Af- 
ter substituting it to the self-consistency equation we get 
the expression for the gap functions 

A,;(fc) ^ R^^N.ik). (6) 

The elements of the matrix R — R{k) are Ru — {XuSi — l) 
and Rij = \jSj, where 




The source functions Nj{k) come from the r.h.s of 
Eqs.®. The strict definition of source functions is given 
in the AppendisEl 

The real space asymptotic of the gap functions ([6]) is 
determined by the contributions of the singularities of 
the response function R~^{k) which are poles at the ze- 
ros of the determinant Dii{k) = Det[i?(A:)] and branch 
points at fc = 2i{lnj/vFj- Similarly to Ref. (jT7l ) we as- 
sume the branch cuts to lie along the imaginary axis from 
k = 2iflnj/vFj to k = ioo. To find the asymptotics of 
the gaps Aj(r), we need only to take into account the 
poles of R~^{k) lying in the upper complex half plane 
below all the branch cuts. In this case all the zeros of 
the function Df({k) are purely imaginary k = z/i„. Each 
of them can be associated with the particular mass ^„ 
of the composite mode formed by a superposition of gap 
functions in two superconducting bands. The compos- 
ite character of the modes arises in our case because the 
two bands arc directly coupled. The inverse of the mass 
controls the characteristic length scale at which this su- 
perposition of the gap fields varies. Therefore the lightest 
mass determines very-long-distance decay of both Ai and 
A2. The contribution from the branch cut contains all 
the length scales which are smaller than the threshold 
one given by position of the lowest branch point k = iq^p 
where q^p = 2 min(17o2/'yF2, ^^oi/^Fi)- 
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The Eq.® results in the asymptotical expression for 
the gap functions 

A,(r) = / driG,,(r,ri)iV,(ri). 

JQ 

Here Nj (r) is the Fourier-Bessel image of the source func- 
tion in Eq.® and 



G{r,ri) = A„/io(g„r)/o(g„ri)- 



(8) 
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where Kq and Iq are MacDonald and modified Bessel 
functions. The matrices An determining the contribu- 
tions of the pole terms are 



dk 



i?22 ^Rl2 
— i?21 Rll 



k-=.9„- (9) 



and the branch cut contribution is determined by the the 
jump of the response function 



R- 



R-^{k^is + Q)- Rr^{k = is-Q). (10) 



Under rather general conditions, the response function 
in Eq.® has two poles given by zeros of the determi- 
nant Dj^{k) = which lie below the branch cuts. Thus 
in this case the asymptotical behaviour of the gap func- 
tions is principally different from the single band super- 
conductor, despite the fact they share the same U{1) 
symmetry of the order parameter. The two poles deter- 
mine the two inverse length scales or, equivalently, the 
two masses of composite gap functions fields, which we 
denote as "heavy" = l^H and "light" = f-L 

(i.e. > Mi)- The corresponding composite gap func- 
tion modes are parameterized by the two "mixing angles" 
9l,6h as follows: 




cosOl smWi 
— sin Oh cos 9h 



Ai 

A2 



(11) 



Note that in the two-band GL theory without interband 
impurity scattering terms one has 9l = G r Be- 
low we recover this behavior at elevated temperatures 
without using GL-like expansion, thereby verifying pre- 
dictions of phenomenological GL models. However, out- 
side the range of validity of the GL theory we find that 
Ol + Oh- 

Let us now consider in detail an example of the system 
with All = 0.25, A22 = 0.213, ui = n2 = 0.5 and various 
values of the interband coupling Aj — X12 = A21. We 
focus on the two different regimes, determined by the 
band parameter "fp = vp2/vFi namely (i) 7i? > 1 and 
(ii) 7F < 1. 

(i) Some of the basic properties of this regime are cap- 
tured by the particular case when = 1. The examples 



of the temperature dependencies of the masses ^l,h{T) 
are shown in the Fig[TlJa). The two massive modes coex- 
ist at the temperature interval < T < Tc, where the 
temperature Tj* is determined by the branch cut position, 
shown in the Figllja) by black dashed line. For temper- 
atures T < Tj* there exists only one massive mode. At 
very low temperatures the mass /^l is very close to the 
branch cut. As the interband coupling parameter is in- 
creased, the temperature Ti rises and becomes equal to 
Tc at some critical value of Aj = Xjc- For the particular 
case of 7i? = 1 we found an exact condition Ajc = A22- 
The evolution of the masses fiL.H is shown in the se- 
quence of plots Figl2lja)-(d) for Aj increasing from the 
small values Aj ^ Aii,A22 to the values comparable to 
intraband coupling Aj ^ Aii,A22- 

(ii) In the case if jp < jth (where jth is a charac- 
teristic value determined by the system parameters) the 
two massive modes coexist at some temperature inter- 
val T* < T < T* where T* < T^. For the particular 
case when jp ~ 0.5, the temperature dependencies of 
IJ'L,h{T) are shown in the Fig[T]^b). 

In Fig[TJa,b) the mixing angles 9l and On given by 
Eg. pT]) are shown by blue dashed and dash-dotted lines 
correspondingly. In the case (i) near the critical tempera- 
ture the angles arc approximately equal, which provides 
for this regime a microscopic verification for of the re- 
sults obtained using phenomenological GL theoriesiSiii. 
At lower temperatures the discrepancy is considerable 
and grows with the increasing interband coupling. Large 
deviations of the mixing angle from and 7r/2 signal 
strong mixing of the gap fields. It occurs near the avoided 
crossing points of ijll{T) and finiT). In case (i) shown 
in FiglIJa) there is one avoided crossing point and in 
the case (ii) in there can be two of them, as shown in 

FigHTb). 

The discussed above existence of two modes associated 
with mixed gap functions can, under certain conditions, 
result in the type- 1.5 behavior as it was demonstrated 
in the framework of GL approacbiSili. However, impor- 
tantly the microscopic formalism we use here allows to 
describe type- 1.5 superconductivity beyond the validity 
of GL models. The type-1.5 behavior requires a den- 
sity mode with low mass to mediate intervortex at- 
traction at large separations, which should coexist with 
short-range repulsion. 

We find that the temperature dependence of /il(T) 
is characterized by an anomalous behavior, which is in 
strong contrast to temperature dependence of the mass 
of the gap mode in single-band theories. As shown on 
Figlljc) the function ^,l{T) is non-monotonic with the 
minimum at the temperature Tmin- The minimum is 
close to the crossover temperature where the second su- 
perconducting band becomes active. The maximum is 
located at the temperature Tmm < Tmax < Tc. 

The structure of the composite gap function mode 
shown in the Figjljc) A^ is characterized by the mix- 
ing angle 9l given by Eq. (fTT|) . At the temperature inter- 
val T < Tmax the mixing angle is 9l ~ tt/2. Therefore 
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FIG. 1: Masses fiL,H of the composite gap function fields for 
(a) 7F = 1 and Aj = 0.005, (b) jf = 0.5 and Xj = 0.0025. 
The position of branch cut is shown by black dashed line. 
The mixing angles 6l,h axe shown by blue dashed and dash- 
dotted lines correspondingly, (c) Temperature dependence of 
the mass (black thich curves) and the corresponding 

mixing angle 9l determined by Eq. (|11|) (red thin curves) for 
7f = 1; 2; 5 (solid, dashed and dash-dotted curves). The 
coupling parameters are An = 0.25, A22 = 0.213 and \j — 
0.005. (d) Temperature dependence of the mass ij.l{T) for 
different values of coupling constant A22. 



in this temperature regime, the mode with lightest mass 
consists primarily of the weak band gap A2 (r) with a tiny 
admixture of Ai(r). Note that in this regime the over- 
all behavior of |Ai|(r) outside the long-range asymptotic 
tail has relatively weak dependence on interband cou- 
pling (i.e. at larger distances from the core it has slowly 
recovering tail associated with only tiny suppression rel- 
ative to its ground state value). At the same time the 
recovery of | A2 1 (r) to a larger degree is dominated by the 
light mass mode. 



C. The high temperature limit. 

As noted above at elevated temperatures the mixing 
angles have close values, consistently with the type-1.5 
behaviour which appears in the framework of two-band 
Ginzburg-Landau models^i. At very high temperatures 
Tmax << T < Tc the mixing angle 9l gradually becomes 
small 0L <C TT, which means that there the mode A^ is 
dominated by the strong band contribution Ai. 

Since any Josephson interband coupling breaks the 
symmetry of the system in question down to U{1), then 
according to Ginzburg-Landau argument this symmetry 



FIG. 2: Masses fiL and fiH (red solid lines) of the compos- 
ite gap function fields for the different values of interband 
Josephson coupling Aj and 7f = 1. In the sequence of plots 
(a)-(d) the transformation of masses is shown for Aj increas- 
ing from the small values Aj <^ Aii,A22 to the values com- 
parable to intraband coupling Aj ~ An, A22. The particular 
values of coupling constants are An — 0.25, A22 ~ 0.213 and 
Aj = 0.0005; 0.0025; 0.025; A22 for plots (a-d) correspond- 
ingly. The branch cuts are shown by black dashed lines. In 
(a) with blue dash-dotted lines the masses of modes are shown 
for the case of Aj = 0. Note that at Aj = the two masses 
go to zero at two different temperatures. Because 1/ ^l,h are 
related to the coherence lengths, this reflects the fact that 
for C/(l) X (7(1) theory there are two independently diverging 
coherence lengths. Note that for finite values of interband 
coupling only one mass goes to zero at one Tc. 



dictates that, asymptotically, in the limit T Tc one 
should recover a single-component-like GL temperature 
dependence fiL ^ \/^ — T/Tc of a single order parameter 
(at the level of mean-field theory }i. 

In the regimes corresponding to Fig[TJa,c) very close to 
Tc the mixing angle of the heavy mode is small 9h ^ 
which makes the contribution of the smaller gap A2 to 
the heavy mode the dominating one. This behaviour of 
the mixing angles, and the fact that that for non-zero 
Josephson coupling only one mass fJ.L{T) goes to zero at 
T ^ Tc allows one to neglect the heavy mode and con- 
struct a mean-field GL order parameter with the scaling 
fiL ^ -^/l — T/Tc as an "asymptotic" characteristic in 
the limit T ^ Tc- However as shown in the Fig[ljc) 
the temperature region of such behavior shrinks drasti- 
cally for large disparities of the band characteristics and 
weak interband couplings. In general the smaller is the 
interband coupling, the closer to Tc one should be in 
order to obtain single component-like GL scaling. For 
a wide range of parameters the mean field GL theory 
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with the single component-like scaling ^il ^ -^/l — T/T^ 
will emerge only inf initesimal ly close to Tc.Note that the 
limit where /il ~ -^/l — T/T^ is in certain cases unphysi- 
cal because the underlying mean-field theory can become 
invalid because of fluctuations, at temperatures lower 
that the temperature where this scaling would take place. 
Thus even in weak-coupling two-band systems with U{1) 
symmetry, for a wide parameter range, one could not ap- 
ply a leading order in (1 — T/Tc) GL theory since the re- 
gion of its applicability will fall into the parameter space 
where underlying mean field theory is not valid because of 
fiuctuations. In contrast to single-component systems, as 
the consequence of the presence of two gaps even slightly 
away from Tc the behaviour of ^J-hiT) can be drastically 
different from the usual GL scaling. As a result the prod- 
uct A/ii where A is the magnetic field penetration length 
acquires a strong temperature dependence. Moreover as 
we show below, its limiting value at Tc does not deter- 
mine entirely the intervortex interaction potential nor the 
magnetic response of the system. Therefore one cannot in 
general parameterize the magnetic response of two-band 
systems by the single GL parameter k = A/^. 



D. Light mode of gap function field and type-1.5 
beliavior. 

The plots of IJ,l{T) for — 1;2;5 are shown in 
FigHJc) by solid, dashed and dash-dotted thick black 
lines. There is a clear general tendency of increasing 
Tmax with growing parameter 7^? which characterizes 
band disparity. It leads to broadening of the tempera- 
ture region of the anomalous behavior of the mass Hl(T) 
where the fields asymptotics are dominated by the weak 
band. The Figlljc) clearly demonstrates the considerable 
overall suppression of fi^ with growing parameter jp. 
The inverse of the mass of the light composite gap mode 
fiL sets the range of the attractive density-density con- 
tribution to intervortex interaction. Therefore the condi- 
tion for the occurrence of the intervortex attraction will 
be met if is smaller than A~^. 

Thus a physically important situation arising in a two- 
band superconductor, is that for a wide range of param- 
eters even slightly away from Tc the temperature depen- 
dence of fiL, is dramatically different from that of the 
inverse magnetic field penetration length A~^. 

Furthermore because the softest mode with the mass 
/ii in two band system may be associated with only a 
fraction of the total condensate, and because there could 
be the second mixed gap mode which can have larger 
mass /.iH, the short-range intervortex interaction can be 
repulsive. Since ultimately the sign of the long range in- 
teraction is decided by the competition of A^^ and /il 
we plot their temperature dependencies in Figl31[a). It 
shows how in these cases the system goes from type-II 
to type-1.5 behavior as temperature is decreased. The 
type-1.5 behavior sets in when /Zi becomes smaller than 
A"-'^, and, the density associated with the light mode is 



small enough that the system has a short-range intervor- 
tex repulsion. 

To contrast the physics of fundamental modes in two- 
band case with singe-band case we plot on FiglSjb) the 
product of A and /i^ . Note that only infinitesimally close 
to Tc, this product can be interpreted as GL parame- 
ter K because the inverse mass ^/2fj,J^^ becomes the sin- 
gle component-like GL coherence length. However away 
from Tc it represents a mass of the softest of competing 
modes and the product Afj.^ has a strong and nonmono- 
tonic temperature dependence shown on Fig. ^h) . 



III. SELF-CONSISTENT CALCULATION OF 

THE VORTEX STRUCTURE AND 
NON-MONOTONIC VORTEX INTERACTION 
ENERGY 

Next we calculate self-consistently the structure of iso- 
lated vortex for different values of 7^- In these calcula- 
tions wc fix the values of parameters ai by adjusting the 
partial DOS which in the case of cylindrical Fermi sur- 
faces is regulated by the ratio of effective masses so that 
n2 = ni/jp and A12 = ^21/ If- We chose the follow- 
ing values of the coupling parameters An = 0.25, X22 = 

0. 213. The interband interaction is small A21 = 0.0025 
and the temperature is T = 0.6 when Aio 3> A20. In 
this case the composite gap function mode AL(r) con- 
sists mainly of the weak gap A2(r). Thus, although at 
the very long ranges the behavior of both |Ai|(r) and 
|A2|(r) are determined by the same mass /i^, the over- 
all behavior (i.e. outside asymptotic regimes) of the gap 
|Ai|(r) [shown by red dashed lines in Figl^l^c)] is not very 
sensitive to the parameter jp. A complex aspect of the 
vortex structure in two-band system is that in general the 
exponential law of the asymptotic behavior of the gaps 
is not directly related to the "core size" at which gaps 
recover most of their ground state values. We can char- 
acterize this effect by defining a "healing" length Lai of 
the gap function as follows \Ai\{LAi) = 0.95Aio. Then 
we obtain that Lai ~ 0.8 for all values of 7^- On the con- 
trary, the healing length La2 of changes significantly such 
that La2 = 1-6; 2.5; 3.2; 3.9; 4.5 for 7^ = 1; 2; 3; 4; 5 
correspondingly. 

To demonstrate the type-1.5 behavior wc have chosen 
the parameters in the self-consistency equation for the 
current such that the characteristic magnetic field local- 
ization length Lh ~ 2 is much larger than Lai- This 
leads to a existence of regular vortex lattices in a wide 
range of strong magnetic fields (i.e. when vortices are 
closely packed and thus experience only strong short- 
range repulsive interaction). However, the high magnetic 
field behavior notwithstanding, the vortex structures 
shown in Figl3{c) clearly shows that Lai Lh <C La2 

1. e. the long-range interaction is attractive and thus the 
system in fact belongs to the type-1.5 regime. 

Next, to demonstrate the type-1.5 superconductivity 
i.e. large-scale attraction and small-scale repulsion of 
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(a) (b) 




FIG. 3: (a) Masses and fin (red solid and dotted lines) 
of the composite gap function fields and inverse London 
penetration (blue dashed lines) for the different values of 
Ajj.L{Tc) /V2 = 1; 2; 3; 5. The position of branch cut is shown 
by black dash-dotted line, (b) The temperature dependence 
of the quantity Afj,L for \^l(Tc)/\/2 — 1;2;3;5 (red solid, 
blue dashed and black dash-dotted lines), (c) Distributions 
of magnetic field H{r)/H(r = 0), gap functions |Ai|(r)/Aio 
(dashed lines) and |A2|(r)/A2o (solid lines) in a single vor- 
tex for the coupling parameters An = 0.25, A22 = 0.213 
and A21 = 0.0025 and different values of the band parameter 
7f = 1;2;3;4;5. (d) The energy of interaction between two 
vortices normalized to the single vortex energy as function of 
the intervortex distance d. In panels (c,d) the temperature is 

r = o.6. 



vortices which originates from disparity of the variations 
of two gaps, we explicitly calculate the intervortex in- 
teraction energy. We evaluate the two-band generaliza- 
tion of the Eilenberger expression for the free energy of 
the two vortices positioned at the points tr = (c?/2,0) 
and tl = (— d/2,0) in xy plane. Here we generalize 
to two-band theory the method developed for calcula- 
tion of asymptotic vortex interaction in singe-component 
theoriesi^. The method assumes that for large separa- 
tions, in the region a; < the fields H, Ai_2 and f[^2 
correspond to the single vortex placed at the point tl 
weakly perturbed by the presence of the second vortex. 
The interaction energy can be expressed through the in- 
tegral over the line x = Q passing in the middle between 



vortices E^t = 2 /_!^ dyEint (y) where 

/oo 
dyH,,Q,+ (12) 
'OO 

The detailed derivation of the above expression can be 
found in the Appendi:xlBl The indices R{L) correspond 
to the solutions of Eilenberger Eqs.(IT]) for isolated vor- 
tices positioned at the points Tj^^l). The first term in 
the Eq.(fT2|) contains the magnetic field Hy{\r — rz,\) and 
the axial component of superfluid velocity distribution 
(5u(|r — tlI) corresponding to the isolated vortex placed 
at the point r = tl. 

In FiglSjd) the interaction energy Eint is shown as a 
function of the distance between two vortices d. The en- 
ergy Eint is normalized to the single vortex energy Ey. 
The plots on Figl31[d) clearly demonstrate the emergence 
of type-1.5 behavior when the parameter is increased. 
This is manifested in the appearance non-monotonic be- 
haviour of Eint{d). 

IV. LOW TEMPERATURE VORTEX 
ASYMPTOTICS AND INTRINSIC PROXIMITY 
EFFECT. 

Finally we discuss the two-band superconductor with 
A20 ^ Aio at T — >■ 0. Note that qualitatively sim- 
ilar regime is realized in the two-band superconductor 
MgB-J^. To model such situation we choose the cou- 
pling constants An = 0.25, A12 = A21 = Aj = 0.05 and 
consider various values of A22. The temperature depen- 
dencies of the mass ij.l{T) for different values of A22 are 
shown in the Fig[ljd). Note that in this case, decreas- 
ing of intraband coupling A22 leads to the decreasing of 
the 11 L Sit low temperatures. This anomalous behaviour 
of the characteristic length scale is clearly manifested in 
the vortex structure shown in Fig.(|4]). The near-core 
gap function profiles [Fig|3Ja,c)] feature shrinkage of the 
vortex core at decreasing temperature, similarly to clean 
single-band superconductors^^. However the asymptotics 
of gap functions [Fig|31[b,d)] are drastically different from 
the single-band case. Indeed, it can be seen that in a 
certain temperature domain the lower the temperature, 
the slower is the recovery the gap functions at large dis- 
tances from the core. Such behavior in the two-band 
system is clearly in a sharp contrast with the overall vor- 
tex core shrinking with decreasing temperature in clean 
single-band superconductors. 

Note that in the above case, at low temperatures we 
have fiL ~ 2^^/K^q + {ttT)^ /vF2- For the especially in- 
teresting regime of purely interband proximity effect- 
induced superconductivity in the weak band we can con- 
sider the limit T ^ A2o/7r. Then « Cjv^i where 
— ^^F2/(27rT) is the coherence length in a pure nor- 
mal metal^ describing the penetration length of super- 
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FIG. 4: Gap function profiles around vortex core for An — 
0.25, A22 = 0.1, A12 — A21 = 0.05. (a) and (c): Variation of 
the gap functions |Aj|(r)/Ajo (j = 1,2) near the core, (b) 
and (d): The behaviour of gap function deviations from the 
vacuum state SAj{r) = 1 — |Aj|(r)/Ajo at longer range. Note 
that in this temperature span the higher is the temperature 
the faster is the long distance decay of SAj{r), which reflects 
the found fact in two-band system the field mass can ibcrease 
with raising temperature [see also Fig[ljd).] 



conducting correlations induced by the proximity effect in 
superconductor/normal metal (SN) hybrid structures'^. 
Thus we obtain that the intrinsic proximity effect due 
to the intcrband couplingii can in certain cases be de- 
scribed by the similar length scale as the usual one 
in SN hybrid structures. At the temperature interval 
A20 <C ttT <C Aio the mass Hl{T) grows linearly with 
temperature [FiglTJd)]. 



V. CONCLUSION 

In conclusion, the rapidly growing family of discov- 
ered multiband superconductors {MgB2, Iron pnictidcs 
etc) requires understanding and classification of possi- 
ble magnetic response of systems with multiple super- 
conducting gaps. Here we reported a microscopic theory 
of magnetic response of a superconductor with two bands 
(the developed approach can be generalized to the case 
of a higher number of bands) . We have shown that new 
physics which arises in multiband systems is the exis- 
tence of several mixed gaps modes. This, in a range of 
parameters results in the existence of the type- 1.5 su- 
perconducting regime. We described the system prop- 



erties and emergence of type- 1.5 regimes in the entire 
temperature regimes, in particular beyond the validity 
of a two-component GL theory. The universal feature of 
all the regimes supporting type-1.5 behavior is the ther- 
modynamic stability of vortex excitations in spite of the 
existence of a mode which varies at a fundamental length 
scale larger than the magnetic field penetration length. It 
results in non-monotonic vortex interaction and appear- 
ance of the additional Semi-Meissner phase in low mag- 
netic fields which is a macroscopic phase separation into 
(i) domains of two-component vortex state and (ii) vortex 
clusters where one of the components is suppressed. 
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Appendix A: Asymptotical behaviour of the gap 
functions. 

Wc focus on the structure of the isolated axially sym- 
metric vortex in two-band superconductor characterized 
by the non-trivial phase winding of the gap functions: 

Ai,2 = |Ai,2|(r)e*'^. (Al) 

We begin by considering the asymptotical behaviour of 
the gap functions at distances far from the vortex core 
when the deviations of all fields from the homogeneous 
values are small. In this case the Eilenderger Eqs.(IT]) can 
be linearized in order to find the asymptotical behavior 
of the gap functions modules |Ai^2|('')- To compare 
with the different linearlization problem in single-band 
case see Ref.©. 

To determine the asymptotic behaviour we use the 
transformation / —J- /e"^, — Z+e"*"^ and rewrite the 
Eilenberger Eqs.(IT]) in terms of the deviations from the 
vacuum state values Aj = A^o — |Aj| and fj = fjQ — fj, 
= — fj~. Then keeping the first order terms 
/s(d) ^iid Aj in the l.h.s. we can rewrite the Eilenberger 
Eqs. in the following form (we omit the band index for 
brevity) : 

vpnVh + 2oj„fd - As (A2) 
z;^nV/, + 2^/s - z^nQ - = X,. 

U)yi ^^n ^^n 

where the higher order terms in Aj, f and are incor- 
porated in the r.h.s. functions ^^(c;) — As(c;)(np,a;„,r). 
In Eas. (jA2[ ) we introduce r2„ ~ y^uj^ + Ag and the func- 
tions /s = / + /+ and fd = f- /+. The higher 
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order terms are incorporated in the functions ^s(d) = 
-'^E(d)(np,t^„,r). 

Then we take the real part of the Eqs.(jA2]) to obtain 
the following system 



(A3) 



4w„ 



Here omit the band index for brevity and denote /^(d) ~ 

Rsf 12(d)- Below we will find the asymptotic of the 
gap fields treating the nonlinear terms in the r.h.s. of 
Eqs. ljASP as source functions. 

The solution of Eas. (jA3[) can be found in the momen- 
tum representation /£ ^(k) = / /£ ^(r) exp(— ikr)d^r. 
Then we get 



8A 



n„ Aftl + (vFk)2 



M(vFk, w„) 



(A4) 



where the last term incorporates the higher order correc- 
tions: 



A/(vFk,w„) = 



Ani + (vpk) 



After substituting it to the self-consistency Eq.([2]) we 
get the expression for the order parameter 



A,;(fc)=i?-.liV,(fc) 



where 



iV,(fc) 



E 

n=0 



Mjd0p 



(A6) 



(A7) 



and the elements of the matrix R — R{k) are defined by 
Ru = [KiS^ - 1) and i?y = >^ijSj, where 



, ,2 /.27r 

^ =4T V ^ / 



(AS) 



n=0 - "J •'U "-ni 

The integrals entering the expressions (|A8[) above are 
'■2'^ dB^ 27r 



62 + (sin6lp)2 h^fWT\ 



so that 



5,(fc)=47rT^ 



y^4172, -I- («i.,fc)2 



(A9) 



The source functions Nj(k) come from the nonlinear 
terms X£ ^ in Eilenberger Eas. (|A3[) . 

The Ea. (|A6p is the two-band response function. To 
compare with the single-band response function see^i. In 
general the real space asymptotic behaviour of the order 
parameter (|A6|) is determined by the contributions of the 



singularities of the response function R~-^{k) which are 
poles and branch points at fc = 2ifl„j/vpj. Analogously 
to the consideration in Ref. (fl7l ) we assume the branch 
cuts to lie along the imaginary axis from k = 2ifl„j/vFj 
to k = ioo. The poles are determined by the zeros of the 
determinant DTi{k) = Deti?(fc) = 0, so that 

DR{k) = (1 - AllSl)(l - A22^2) - Ai2A2l5l^2. 

Since we are interested in the asymptotic behaviour of 
the order parameter, we need only to take into account 
the poles of R~^{k) lying in the upper complex half plane 
below all the branch cuts. In this case all the zeros of the 
function Dji{k) are purely imaginary k* = ig„. Each 
of them can be associated with the particular mass of 
the gap function field ~ ^/Qu which determine the 
characteristic length scale of the gap function variation. 
On the other hand the contribution from the branch cut 
contains all the length scales which arc larger than the 
threshold one given by position of the lowest branch point 
k — iqhp where 



Qbp = 2min(Oo2/wF2,^^oi/wfi)- 



(AlO) 



Appendix B: Energy of interaction between two 
vortices 

1. General free energy expression 

The two-band generalization of the Eilenberger expres- 
sion for the free energy^ reads as follows 



F(r) 



-pillAip -^-p22|A2| 



(Bl) 



where 



pj (AiA* + A2A^) + Fn+Fj2 



Pll Pi2 \ ^ f Pll Pl2 

P21 P22 J \P21 P22 



P.J = Pl2 = P21 



and 



Fjj = 2 y2 / njlj{ujn,0p,r)d9p 



with 



/,(w„,0p,r) = A*/,+A,/; 



(B2) 



-I- [gj - 1) 2aj,; 
where j = 1,2 and 



^npV (In/, -In/; 



UJn = LOn + iVpjXlpA/2. 

Then the variation of the free energy (|Bip with respect to 
the fields A and A gives the self-consistency Eqs.Q and 
(PI correspondingly. The variation over / and /+ with 
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the normalization condition taken into account yields the 
Eilenberger Eqs.([l]). Provided the functions f,f~^,g sat- 
isfy the Eqs.([T]) the expression (jB2| can be rewritten as 



/j(w„,6'p,r) 



(B3) 



The first term in this expression cancels with the second 
and forth terms in Ea. (|B5p . For the variation of magnetic 
field energy in Eq. (|Bip we obtain 

Hvi5H = V • ((5Q X Hv) + V X Hv • <5Q. 

Then we are left with the non-zero terms 



2. Linearized theory of vortex interaction 

To calculate the energy of vortex interaction we eval- 
uate the free energy expression for the system of two 
vortices positioned at the points tr = (d/2, 0) and 
tl = (— d/2,0) in xy plane. Here we employ the method 
similar to that ini^. 

Let us consider the half-plane a; < containing only 
one of the vortices. We decompose the gap function 
into amplitude and phase (we omit the band index for 
brevity) 



A(r) = |A|(r)exp(i$). 



(B4) 



The total phase can be written in the following form $ 



^ns, where 



*L(H)(r) = arctan 



/ y - VRjL) \ 



are the vortex phases and $,is(r) is a regular part of 
the phase. At the region x < we can make the gauge 
transformation removing the phase $fl(r), since it does 
not contain singularities. After this transformation we 

can assume that the fields A, Ai^2 and f[~^^ correspond 
to the solutions for a single vortex placed at the point tl 
weakly perturbed by the presence of the second vortex. 



A = Av + SQ; Aj = A 



V] 



where we have introduced the superfluid velocity induced 
by the second vortex SQ = Ar — W^b- Then we obtain 



6L 



{5A,f+ + 5A*U,) + (A,„,5/+ + A://,)(B5) 
+ivFj{gvj - l)np(5Q + 2a)„(55j 
+VFj'-fn^V (ln/„,-ln/+) 



-VFj- 



2 "P^ I /„j 



7^ 



where 



Sg,^-{f,,Sf+ + f+5f,)/2g,,. 

The last two terms in Eq. (|B5p can be rewritten as follows 

1 
2g^ 
("pV) 



[<S/(npV)/+-5/+(npV)/„] 



(B6) 



{9v - 1) 



Sf Sf^ 



fv fv 



(5F = V • (5Q X Hv 



2^ E ''i^'Fi 



2tt 




The energy of vortex interaction is Eint = 2 / SFdr. It 
can be expressed through the integral over the line x = 
so that Eint = 2 dyx ■ Si^t 



e-int = X Hv- 



(B8) 



dOpn^ 



{gv 




To evaluate the second term in Eq. (jB8|) it is convenient 
to bring the Eqs.(IT]) to the gauge invariant form3^ decom- 
posing the gap functions into amplitude and phase (jB4[) 
and transforming the Green's functions as / — >■ /e**, 
/+ /+e^'*. Then at the line a; = we can put 



Also we denote Sfj 



where foj — ^uj/y^oj v ^n- ^i^w wv. u^^wui. ujj 

fuj, Sfj' = f^j [L{R) stand for left (right) vortices] 
Therefore up to the second order terms we obtain 



i9vj - 1) 



^ fvj fvj / 



- /A,) (B9) 



2.90 J 

go, - 1 

J 03 

Now we use the symmetry relations fL.Rii^xTn-y) ~ 
f*R,L{-'<^x,ny) and /*(-np) = /+(np). Then the con- 
tribution to the interaction energy (jB8[) from the first 
order term in Eq. ljBOp cancels with the analogous con- 
tribution from the left vortex. Also from the symmetry 
relations we obtain 

Re / cos0pfLfRd0p^O (BIO) 
Jo 

Re cosOpf+f+dOp^Q. 







On the other hand 



/>Z7r 

Im / cos 9p {fLfR - fUi) dOp = 0. 
Jo 
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Therefore we get for the interaction energy and Uj = k '^rijVpj. 

/oo 
-oo 

where 

Eint = HyQy+ (Bll) 

7 = 1,2 w„>0 " 
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